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Abstract 

The model-independent determination of the Higgs-boson CP properties at the Photon 
Collider at TESLA has been studied in detail, for masses between 200 and 350 GeV, using 
realistic luminosity spectra and detector simulation. We consider a generic model with 
the CP violating Higgs tensor couplings to gauge bosons. We introduce a new variable 
describing angular distributions of the secondary WW and ZZ decay products which is 
very sensitive to the CP properties of the Higgs-boson. Understanding of the detector 
performance turns out to be crucial, as the influence of the acceptance corrections is 
similar to the effects of CP violation. From the combined measurement of invariant mass 
distributions and various angular distributions the phase describing a CP violation can 
be determined to about 50 mrad after one year of Photon Collider running. 

1 Introduction 

The physics potential of a Photon Collider [1] is very rich and complementary to the physics 
program of the e + e~ and hadron-hadron colliders. It is an ideal place to study the mechanism 
of the electroweak symmetry breaking (EWSB) and the properties of the Higgs-boson. 

In paper [2] we performed a realistic simulation of the Standard Model (SM) Higgs-boson 
production at the Photon Collider for W + W~ and ZZ decay channels, for Higgs-boson masses 
above 150 GeV. From the combined analysis of W + W~ and ZZ invariant mass distributions 
the 77 partial width of the Higgs boson, T 77 , can be measured with an accuracy of 3 to 8% 
and the phase of 77 — > h amplitude, <^ 77 , with an accuracy between 30 and 100 mrad. In paper 
[3] we extended this analysis to the generalised Standard Model-like scenario of the Two 
Higgs Doublet Model II, 2HDM(II), with and without CP- violation. In case of the solution 



B h with a weak CP violation via H — A mixing, the mixing angle $ha can be constrained 
to about 100 mrad for low values of tan/3. We found that for a general 2HDM (II) with CP 
violation, only the combined analysis of LHC, LC and Photon Collider measurements allows 
for the precise determination of Higgs-boson couplings and of CP- violating H — A mixing angle 

[4]- 

In this paper our aim is to establish CP-property of the Higgs-boson in a generic model with 
a direct CP-violation. The measurement of the Higgs-boson properties at the Photon Collider 
at TESLA is studied in detail for masses between 200 and 350 GeV, using realistic luminosity 
spectra and detector simulation. The model with generic, CP-violating Higgs-boson couplings 
to vector bosons [5] leads to different angular distributions for a scalar- and pseudoscalar-type of 
couplings. From a combined analysis of invariant mass distributions and angular distributions 
of the W + W~ and ZZ decay-products the CP-parity of the observed Higgs state can be de- 
termined independently on a production mechanism [6]. Results given in this paper supersede 
results presented in the second part of our earlier work [7]. 



2 Generic Higgs model with CP violating couplings 

Following the analysis described in [5] we consider a generic model with a direct CP violation, 
with tensor couplings of a Higgs boson, H, to ZZ and W + W~ given by: 

9hzz - ig ^— I Xh • g + Xa ■ £ > 

9hww = i 9 M w (A* ■ <r + A, • sT" [Pl+P '\^- n) ° ) , (1) 

where pi and P2 are the 4-momenta of the vector bosons. Contributions to gauge couplings 
with X H have a structure of the CP-even SM Higgs boson coupling, 1 whereas the ones with 
correspond to a general CP-odd coupling for the spin-0 boson. Coefficients Xh and A^ can be 
written as: 

X H = A • cos $cp , 

A A = A • sin$ CP . (2) 

The couplings of the Standard Model Higgs boson are reproduced for A = 1 and §cp — 
(i.e. \h = 1 and Xa = 0). Below we will limit ourselves to A ~ 1 and \$cp\ ^ 1 region, 
corresponding to small deviation from the Standard Model coupling. However, we do not make 
any assumptions concerning Higgs-boson couplings to fermions and we allow for deviations from 
SM predictions in T 77 and 77 . Therefore our results do not depend on the assumed Higgs- 
boson production mechanism and our approach can be considered as model-independent. To 
simplify the analysis, we only assume that the 7i branching ratios to W + W~ and ZZ are the 
same as in the Standard Model. No substantial deviations in the W + W~ and ZZ branching 
ratios are expected unless Higgs-boson Yukawa couplings to up- or down-type fermions become 
very large. Possible small deviation in the branching ratios are equivalent to T 77 variation. 

1 Othcr possible CP-even tensor structure, (pi +P2) M (pi give the angular distributions similar to that 

of the SM Higgs boson and therefore we will not consider this case separately. 
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3 Angular distributions for secondary decay products 



Here we consider the Higgs boson decays to gauge bosons. The angular distributions of the 
secondary W + W~ and ZZ decay products turn out to be very sensitive to the CP properties of 
the Higgs-boson [5]. Angular variables which can be used in the analysis are defined in Fig. 1. 
To test CP-properties of the Higgs-bosons the distributions of the polar angles Gi and G 2 as 
well as the A0 distribution, where A<j) is the angle between two Z- or two W-decay planes, are 
used. Here we propose to consider, instead of two-dimensional distribution in (cosGi, cosG 2 ) 
the distribution in a new variable, defined as 



The (^-variable corresponds to the ratio of the angular distributions expected for the decay 
of a scalar and a pseudoscalar (in a limit y> Mz) [5]. It proves to be very useful and 
complementary to A</> variable. 

The angular distributions in A(p and £, expected for a scalar ($c*p — 0) and a pseudoscalar 
($ CP = |) higgs decays, Ti — > ZZ — > l + l~jj, are compared in Fig. 2. They clearly distinguish 
between decays of scalar (if) and pseudoscalar (A) higgs. For the scalar higgs the distributions 
are almost flat in both Acp and (. For the pseudoscalar case these two distributions have totally 
different character. From the measurement of one or other, or both of these distributions it's 
possible to establish the CP properties of the Higgs boson, even without taking into account 
the production mechanism. 

Note, that the theoretical distribution of the higgs decay-angle G-^ is expected to be flat 
both for scalar and pseudoscalar higgses, as this is characteristic to spin particle. However, 
the measured G^ distributions for if and A decays can differ due to the different acceptance 
corrections and interference with a non-resonant background. Therefore, in the observed G^ 
distributions, as well as reconstructed 77 invariant mass distributions for W + W~~ and ZZ 
decays, some sensitivity to the Higgs-boson couplings can show up. 



4 Reconstruction of the angular distributions 

Our analysis uses the CompAZ parametrisation [8] of the realistic luminosity spectra for a 
Photon Collider at TESLA. We assume that the centre-of-mass energy of colliding electron 
beams, y/s^, is optimised for the production of a Higgs boson with a given mass. All results 
presented in this paper were obtained for an integrated luminosity corresponding to one year of 
the photon collider running. The total photon-photon luminosity increases from about 600 fb _1 



for yfs^ = 305 GeV (optimal beam energy choice for M n = 200 GeV) to about 1000 fb 1 for 
v /5^ = 500 GeV (M n = 350 GeV). 




Measurement of the angular distributions has been studied using the samples of ZZ events 
generated with the PYTHIA [9] and SIMDET [10] programs, as described in [2]. The resolutions 
of the reconstructed decay angles Q and <fi, for the scalar Higgs-boson with mass of 300 GeV 



c = 



sin 2 Gi • sin 2 Q 2 



(3) 



(1 + cos 2 Gi) • (1 + cos 2 G 2 ) ' 
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Figure 1: The definition of the polar angles Q n , 0i and 2 , and the azimuthal angles 0i and 
(f) 2 for the process 77 — > h — > ZZ — > 4 /. A0 is the angle between two Z decay planes, 
A0 = 02 — 0i- All polar angles are calculated in the rest frame of the decaying particle. 




Figure 2: Normalised angular distributions in A0 (left plot) and £ (right plot), expected for 
scalar and pseudoscalar higgs decays H,A^ ZZ — > l + l~jj, for the higgs mass of 300 GeV. 
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(primary electron-beam energy of 209 GeV) are compared in Fig. 3 for the leptonic and hadronic 
Z decays. 

A very good resolution is obtained for all considered distributions, nevertheless the measured 
angular distributions are strongly affected by the selection cuts used in the analysis. Because 
of the cut on the lepton and jet angles, applied to preserve a good mass resolution (see [2] for 
more details), we observe a significant loss of the selection efficiency for the events with lepton 
or jet emitted along the beam direction. Also the Durham jet algorithm used in the event 
reconstruction imposes constraints on the angular separation between leptons and jets. Both 
effects result in a highly non-uniform angular acceptance. The obtained selection efficiencies 
for 77 — > ZZ — > l + l~jj events, as a function of the angle A<p, are presented in Fig. 4. Addi- 
tional cuts on the reconstructed ZZ invariant-mass have been introduced to optimise the signal 
measurement: for Higgs boson mass of 300 GeV the accepted mass range lies between 286 and 
312 GeV. The efficiencies for ZZ events coming from decays of scalar, pseudoscalar higgs and 
non-resonant background are compared. The polar-angle distributions differ among these three 
classes of events, therefore significantly different acceptances are obtained. Understanding of 
this effect is crucial, as it can mimic a pseudoscalar type of distribution. For the first time such 
effect is taken into account in the analysis of the considered process. 

In the following analysis, the reconstructed Acfi values range from to n, since we are not 
able to distinguish between quark and antiquark jet. 



5 Analysis 

For arbitrary values of model parameters A and §cp we calculate the expected angular and 
invariant mass distributions for ZZ and W + W~ events by convoluting the corresponding cross- 
section formula with the analytic photon-energy spectra CompAZ [8]. To take into account 
detector effects, we convolute this further with the function parameterising the invariant-mass 
resolution and the acceptance function containing the angular- and jet-selection cuts. 2 For the 
measurement of the event distributions in A0 and ( we introduce the additional cuts on the 
reconstructed ZZ invariant mass (for ZZ events) or on the reconstructed W + W~ invariant 
mass as well as on the higgs-decay angle 6^ (for W + W~ events). The cuts were optimised for 
the smallest relative error in the signal cross-section measurement. 

The expected precision in the measurements of the A0- and of the (^-distributions, for 
77 — > ZZ — > l + l~jj events and 77 — > W + W~ — * 4j events is illustrated in Figs. 5 and 6, 
respectively. Calculations were performed for the primary electron-beam energy of 152.5 GeV 
and the Higgs-boson mass of 200 GeV. The results are compared with the expectation of the 
generic model with <& C p = (as in SM) and Q C p — f • For better comparison of the shape 
of distributions, results for pseudoscalar Higgs-boson couplings ($cp — f ) are normalised to 
the Standard Model expectations. We see, that even after taking into account beam spectra, 
detector effects, selection cuts and background influence, the differences between shapes of the 
angular distributions for the scalar and pseudoscalar couplings are still significant. Therefore 
we should be able to constrain Higgs-boson couplings from the shape of the distributions, even 

2 To simplify acceptance calculations, we neglect effects of the finite angular resolution. 
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Figure 3: Resolution in the reconstructed Z-decay angles 6 and 0, for the leptonic (6;, <pi) 
and hadronic (Qj, 4>j) final states. Events were simulated with the PYTHIA and SIMDET 
programs, for a primary electron-beam energy of 209 GeV and Standard Model Higgs-boson 
mass of 300 GeV. 



o 
c 

0) 

'o 

<U 0.4 



0.2 



SM higgs 

pseudoscalar 

background 



6 

gen. 



Figure 4: Selection efficiency as a function of the angle A(p between two Z decay planes, for 
ZZ — > l + l~jj events coming from the scalar higgs decays, pseudoscalar higgs decays and non- 
resonant background. Events were simulated with the PYTHIA and SIMDET programs, for a 
primary electron-beam energy of 209 GeV and Higgs-boson mass of 300 GeV. Only events with 
the reconstructed higgs mass between 286 and 312 GeV are accepted. 
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Figure 5: Measurement of the angle A(f)zz between two Z-decay planes (left plot) and of the 
variable (zz calculated from the polar angles of the Z — ► and Z — > jj decays (right 
plot) for ZZ — ► l + l~jj events. Points with error bars indicate the statistical precision of the 
measurement after one year of Photon Collider running at nominal luminosity. The solid (red) 
and dashed (blue) lines correspond to predictions of the model with pure scalar (<3?cp = 0) 
and pseudoscalar ($op = f ) Higgs-boson couplings, respectively. Green line represents the SM 
background of non-resonant ZZ production. Signal and background calculations are performed 
for primary electron-beam energy of 152.5 GeV and the Higgs-boson mass of 200 GeV. 
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Figure 6: Measurement of the angle Acfrww between two W-decay planes (left plot) and of the 
variable (ww calculated from the polar angles of two W — ► jj decays (right plot) for WW — * 4j 
events. Points with error bars indicate the statistical precision of the measurement after one 
year of Photon Collider running at nominal luminosity. The solid (red) and dashed (blue) 
lines correspond to predictions of the model with pure scalar ($cp = 0) and pseudoscalar 
($ CP = |) Higgs-boson couplings, respectively. Green line represents the SM background of 
non-resonant ZZ production. Signal and background calculations are performed for primary 
electron-beam energy of 152.5 GeV and the Higgs-boson mass of 200 GeV. The insets show 
the relative deviations from the Standard Model predictions ($cp — 0) for pseudoscalar Higgs- 
boson couplings. 
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if the overall normalisation related to the Higgs-boson production mechanism is not known. For 
the Standard-Model Higgs-boson decays to ZZ, about 675 Higgs-boson events and 145 non- 
resonant background events are expected after all selection cuts, in the selected mass window 
from 180 to 210 GeV. For the WW channel about 8000 signal events are expected, compared to 
about 170 000 background events. In both cases a statistical precision of the signal measurement 
is similar. However, due to overwhelming background contribution, the WW analysis is much 
more dependent on systematic uncertainties. 

Each of the considered angular distributions discussed above can be fitted with model ex- 
pectations, given in terms of the parameters A and &cp describing Higgs-boson couplings to 
gauge bosons, parameters T 77 and </> 77 describing the production mechanism, and an overall 
normalisation. Statistical errors in the determination of $cp, resulting from fits to the different 
distributions for ZZ events, are compared in Fig. 7. The remaining model parameters, except 
for the normalisation, are fixed to the Standard Model values. Out of the three considered 
angular distributions, parameter £ turns out to be the most sensitive to the angle §cp describ- 
ing the CP-violation in the Higgs-boson couplings. Surprisingly, the smallest error is obtained 
from the fit to the invariant mass distribution. This is because, due to the different selection 
efficiencies for the scalar and pseudoscalar decays and also to interference effects, the relative 
normalisation of the Higgs-boson signal to the non-resonant background depends on the angle 
$<7P. Therefore, we include the invariant mass distributions for W + W~ and ZZ decays in 
the combined fit. When all parameters are allowed to vary in such a fit, the invariant-mass 
distributions constrain mainly the T 77 and 77 . 




Figure 7: Statistical error in the determination of §cp, expected after one year of Photon 
Collider running, as a function of the Higgs-boson mass M-n- Fits were performed to the shape 
of the three individual angular distributions measured for ZZ events, as indicated in the plot, 
and to the reconstructed invariant mass distribution. Results of the simultaneous fits to three 
angular distributions and to all distributions are also shown. All other parameters, except 
for the overall normalisation, are fixed to the Standard Model values. Errors were calculated 
assuming § C p rs 0. 
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6 Results 



We calculate the expected statistical errors on the parameters A and $c*p ; from the combined 
fit to angular distributions measured for ZZ and W + W~ decays, and to the invariant mass 
distributions. Results are shown in Fig. 8. The two photon width of the Higgs boson, T 77 , the 
phase 77 and normalisations of both samples are allowed to vary in the fit, so the results are 
independent to the production mechanism. One observes that for low Higgs-boson masses below 
250 GeV, better constrains are obtained from the measurement of W + W~ events, whereas for 
masses above 300 GeV smaller errors are obtained from the measurement of ZZ events. The 
error on &cp expected from the combined fit is below 50 mrad in the whole considered mass 
range. The corresponding error on A is about 0.05. 



7 Summary 

An opportunity of measuring the Higgs-boson properties at the Photon Collider at TESLA has 
been studied in detail for masses between 200 and 350 GeV, using realistic luminosity spectra 
and detector simulation. We considered measurement of the invariant mass distributions and 
the various angular distributions. A new variable is proposed to describe the polar angle 
distributions of the secondary W + W~ and ZZ-dec&y products. Event reconstruction and 
selection procedure result in acceptance corrections which are highly non-uniform and depend 
on the Higgs boson CP properties. Understanding of these effects is crucial in the analysis of 
angular distributions. From the model-independent study the angle describing a CP violation 
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Figure 8: Statistical error in the determination of $cp (left plot) and A (right plot), expected 
after one year of Photon Collider running, as a function of the Higgs-boson mass M^. Combined 
fits were performed to the considered angular distributions and invariant mass distributions for 
ZZ events and W + W~ events. Results were obtained assuming small deviations from Standard 
Model predictions, i.e. A ~ 1 and $cp « 0. Two photon width of the Higgs boson, r 77 , the 
phase 77 and normalisations of both samples are allowed to vary in the fit. 
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in the generic Higgs-boson couplings to vector bosons can be determined with accuracy of about 
50 mrad. 
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